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Context. The Local and Loop I superbubbles are the closest and best investigated supernova (SN) generated bubbles and serve 
as test laboratories for observations and theories of the interstellar medium. 

Aims. Since the morphology and dynamical evolution of bubbles depend on the ambient density and pressure distributions, a 
^b : realistic modelling of the galactic environment is crucial for a detailed comparison with observations. 

f^*) ' Methods. We have performed 3D high resolution (down to 1.25 pc on a kpc-scale grid) hydrodynamic simulations of the Local 
\^ , Bubble (LB) and the neighbouring Loop I (LI) superbubble in a realistically evolving inhomogeneous background ISM, disturbed 
C^ • already by SN explosions at the Galactic rate for 200 Myr before the LB and LI are generated. The LB is the result of 19 SNe 
occurring in a moving group, which passed through the present day local Hi cavity. 

Results. We can reproduce (i) the Ovi column density in absorption within the LB in agreement with Copernicus and recent 
1 '• FUSE observations, giving Novi < 2 x 10^^ cm~'^ and Nov/ < 7 x lO^'^ cm~^, respectively, (ii) the observed sizes of the Local 
O ' and Loop 1 superbubbles, (iii) the interaction shell between LB and LI, discovered with ROSAT, (iv) constrain the age of the 
,y ■ LB to be 14.5±o;4 Myr, (v) predict the merging of the two bubbles in about 3 Myr, when the interaction shell starts to fragment, 
^ ' (vi) the generation of blobs like the Local Cloud as a consequence of a dynamical instability. 

. . I Conclusions. We find that evolving superbubbles strongly deviate from idealised self-similar solutions due to ambient pressure 
^ . and density gradients, as well as due to turbulent mixing and mass loading. Hence, at later times the hot interior can break 
l^j ' through the surrounding shell, which may also help to explain the puzzling energy "deficit" observed in LMC bubbles. 

j^ Key words. Hydrodynamics - Shock Waves - ISM: general - ISM: bubbles - ISM: structure - ISM: kinematics and dynamics 

1. Introduction Local Bubble as inferred from a reanalysis of Copernicus 

^ , . 111,. TT • 1- data (Shelton & Cox 1994). One severe restriction in all 

Our solar system is embedded m an Hi cavity extending ,, , , , . , • , ;i t n i. i. ^u i^ r cat 

, „„„ • , ^ 1 . , 1 , , r.r.r^ tlic modcls, which consider the LB to be the result of SN 

about 200 pc m the Galactic plane and roughly 600 pc per- , . . , , • • ^ i i • . 

,. , . ,. ^_. , . ,. ,. explosions, viz. the expansion into a homogeneous ambient 

pendicular to it, according to JNai absorption line studies ,. i u j • ii • i i- ^ i 

, , , , /x 11 1 r,r,r,^\ ^ mcdium, has been removed m the simulations presented 

towards background stars (Laliement et al. 2003). l<or a , . „ .t t -n t ^- ^ ^ -l. c j 

. ° ^ / here. As a corollary, the LB evolution has to be performed 

long time it has been known that at least part 01 this so- • ■ ,1 -^i ^i ^ r j.u t t/ti\ u uui Ti<^oArn 

,, 1 T 1 T^ 111 /TT^N 1 ^■rr ■ ■ ■ r joiutly With that of the Loop I (LI) supcrbubblc. ROSAT 

called Local Bubble (LB) shows dmuse emission m soft ticti/-i u j.- i- u ^ v i j 

, r y.^ ^ ri 1 i r,r.r.\ xn • i FSPC observations, revealing a coherent X-ray shadow 

X-rays (ct. McCammon & Sanders 1990). Ever since the , j t i i j ^ ^i. i- /t^ c a i, u i 

.r „ 1 , 1 1 . p. . ,. . ,. -r„r towards LI, led to the suggestion (Egger & Aschenbach 

LB has been the subiect of intensive studies m radio, UV, -.nnr-x ; i ; ; i. t n • • i ^ -l -ii ^i. j 

^^^^. T r ^r 11 n 1 • 1 1 1995) that the LB is m close contact with the expand- 

EUV and soft X-rays, as well as of analytical and numer- . . , , . t i i i i i t i ^i • • ^ i r 

. , , „. , ■' ' 0,1 „„„ Tx • 1 1 mg neighbouring LI superbubble. LI thus originated from 

ical modelling (e.g. Cox & Anderson 1982, Breitschwerdt ,, . i • • .1 o r^ i i- -^u 

. ^ , , \°, „ . 1 „ ^ „„„, X TT 1 multi-supernova explosions m the Sco Cen cluster, with 

& Schmutzler 1994, Smith & Cox 2001). However, due to ,i ^^ ,1 ^^ , r. . • . n •; 1 , -n ■ j- j 

' ' ' the North Polar Spur being part of its shock illuminated 

the wealth of data and the inherent complex structure of , 1 ,1 

T T, 1, 1 1 1 1 . . outer shell. 



It has also been suggested that part of the soft X- 



the LB, all models have severe shortcomings, concerning 

either the dynamical evolution or the model spectra, in . . 1 1 1 r 1 1 ■ ■ ; 1 1 

1 ,, T i- 1 111 1 ii 1 ray emission could be 01 very local origin, generated by 

some cases both. In particular models based on thermal , , . , 1 . , . . . 

, ^. 11 r -1 i 1 ii 1 11^ charge exchange reactions between solar wind ions and he- 
conduction usually lail to reproduce the observed low Ovi ,. T . / V, r^r^r^r^\ 1 • 1 

, ^. 1 T -^ r 1 ^ 1f^1s -9 • ii liospheric (e.g. Cravens 2000) or even atmospheric plasma, 
absorption column density ol ~ 1.6 x 10 "^ cm ^ m the tt • ;i , r i i r ,, 
However, even m the most favourable case oi all emis- 
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missing flux along other lines of sight, in particular per- 
pendicular to the Galactic plane. 

In this letter we describe for the first time a realistic 
evolutionary scenario for the origin and evolution of the 
Local Bubble jointly with the LI superbubble. The crucial 
differences as compared to previous models are: (i) a re- 
alistic background medium, that has been pre-structured 
by previous generations of SN explosions, which is there- 
fore very inhomogeneous, both in density and pressure 
distribution, (ii) taking into account the time sequence 
and locations of SN explosions, which generate the two 
bubbles, calculated according to the mass distribution of 
extinct members of a moving group and the Sco Cen clus- 
ter, respectively, both derived from a Galactic initial mass 
function (IMF) normalised to still existing low mass stars, 
and (iii) to calculate the density and temperature struc- 
ture of the LB and LI in 3D with high resolution in order 
to sample Ovi along many different lines of sight and com- 
pare them to observations. 

This letter is organised as follows: Section 2 describes 
details of the model setup, in Section 3 our results with 
respect to bubble evolution and comparison to observa- 
tions are discussed, and we close with a general discussion 
and conclusions in Section 4. 
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2. Model and Simulations 

We use a 3D parallelised adaptive mesh refinement (AMR) 
hydrocode to track small scale structures down to 1.25 pc, 
where necessary, and we follow the LB and LI evolution 
in an ISM, driven by SNe types la, Ib-|-c and II occur- 
ring at the Galactic rate (taken from Cappellaro et al. 
1999), on a Cartesian grid of < (x,y) < 1 kpc size in 
the Galactic plane, and —10 < 2; < 10 kpc perpendicular 
to it. This is an extension of the SN-driven ISM model 
of Avillez (2000), fully tracking the time-dependent evo- 
lution of the large scale disk and the Galactic fountain. 
The gravitational field of the stellar disk is adopted from 
Kuijken & Gilmore (1989). Further processes include ra- 
diative cooling, assuming an optically thin gas in coUi- 
sional ionisation equilibrium (CIE, i.e. ionisation rate by 
collisions between electrons and ions is balanced by ra- 
diative recombination rate) using the cooling functions of 
Sutherland & Dopita (1993) for T > 10^ K and Dalgarno 
& McCray (1972) for T < lO"* K, with a temperature cut- 
off at 10 K, and uniform heating due to starlight varying 
with height z (Wolfire et al. 1995). Galactic OB associa- 
tions form in regions with density n > 10 cm""^ and tem- 
perature T< 100 K, with the number, masses and main 
sequence life times of the stars in the association being de- 
termined from the IMF. All OB stars are allowed to drift 
away from their parent associations with random veloci- 
ties assigned according to observations, before exploding 
either in clusters or in the field with a canonical energy of 
lO'"^! erg. 

The LB itself is carved by 19 successive explosions of 
the subgroup Bl of Pleiades, of which some low mass stars 
are now Sco Cen cluster members according to the model 
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Fig. 1. Top: Colour coded temperature map in the range 
10 < T < 10^ K for a slice through the data cube (repre- 
senting the Galactic midplane) of a 3D LB high resolution 
simulation representing the present time (i.e. 14.5 Myr af- 
ter the first explosion) with the LB centred at (175, 400) 
pc and LI shifted 200 pc to the right. Bottom: Same, show- 
ing the "future" of the LB and LI at t = 29.7 Myr. Note 
that part of the LB has merged with the ISM, part has 
been engulfed by LI. 



of Berghofer & Breitschwerdt (2002, henceforth BB02). 
The adjacent LI bubble is the result of explosions of Sco 
Cen cluster members (see Egger 1998). Applying an IMF 
for OB associations (Massey et al. 1995) and following 
BB02, we find that about 39 SNe have occurred in LI until 
now, and 38 SN candidates are expected to explode within 
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the next 13 Myr. In contrast, the LB is extinct now, the 
last explosion having occurred about 0.5 Myr ago. Main 
sequence life times are calculated from Stothers (1972). 

As initial conditions for the present run we used data 
cubes of an evolved SN driven ISM, in which global dy- 
namical equilibrium has been established (cf. Avillez & 
Breitschwerdt 2004 (henceforth AB04), 2005a). As the 
ISM structure reveals a typical pattern size of a few 
hundred pc (Avillez & Breitschwerdt 2005b, henceforth 
AB05), we chose arbitrarily a site with a sufficient amount 
of cold HI mass to form the Sco Cen cluster, here located 
at (375,400) pc (cf. Fig.[TIl. Then we followed the trajec- 
tory of the Pleiades moving group Bl, whose SNe in the 
LB went off along a path crossing the location at (175, 400) 
pc, chosen to match the present observed distances from 
the Sun and the Sco Cen. Due to small peculiar motions 
of the local gas, the LB will be tied to the local standard 
of rest (LSR; for details of the moving group kinematics 
see BB02). Periodic boundary conditions are applied along 
the four vertical boundary faces, while outflow boundary 
conditions are imposed at the top (z = 10 kpc) and bot- 
tom (z == —10 kpc) boundaries. The simulation time of 
the LB and LI evolution run was an additional 30 Myr. 

3. Results 

3.1. Evolution of the Local and Loop I Superbubbles 

The present runs were started after global dynamical equi- 
librium, specifically the Galactic Fountain, had been es- 
tablished (i.e. after 200 Myr of evolution for a Galactic SN 
rate; for details, see AB04). The locally enhanced SN rates 
produce coherent LB and LI structures within a highly 
disturbed and still evolving background medium (due to 
ongoing star formation). Successive explosions heat and 
pressurise the cavities, which at first look smooth, but de- 
velop internal structure after some time; for the LB this 
occurs at i > 8 Myr after the first explosion. After ~ 14.5 
Myr (present time) 19 SNe have exploded inside the LB 
cavity, which becomes progressively elongated and reaches 
a size of 180 by 220 pc in the Galactic midplane (Fig. ^ 
top), due to the increasing influence of ambient density 
and pressure gradients. Although magnetic flelds and cos- 
mic rays have an influence on bubble evolution by tension 
and pressure forces, it seems unlikely that our results will 
change significantly, because in the LB the last SN went off 
only 0.5 Myr ago, and Loop I is still active. If the bubble 
pressure decreases markedly there may however be some 
deviations. Both bubbles are bounded by shells whose in- 
teraction already generates Rayleigh- Taylor instabilities 
due to a larger pressure in LI with respect to the LB, in 
agreement with a linear stability analysis of Breitschwerdt 
et al. (2000). After becoming nonlinear this will lead to the 
formation of cloudlets, which will travel towards the solar 
system from the Sco Cen region, in agreement with obser- 
vations of local clouds with measured velocities of about 
26 km/s. On the basis of our simulations we further pre- 
dict that the interaction shell will break up in 3 Myr from 



now, allowing mass transfer of hot gas from LI to the LB, 
and in '^ 15 Myr from now the bubbles will have merged 
(see Fig. n bottom). 

In CIE, Ovi is most abundant at T - 3 x 10^ K, 
whereas soft X-Ray emission of a thermal plasma requires 
a temperature of ^ 10^ K. Fig. |21 shows the temperature 
distribution in the critical range, 10^ < T < 10^ K, dur- 
ing the critical time, 14.5 < t < 15.0 Myr, when the LB 
plasma starts cooling down. It is obvious that between 14.5 
and 14.7 Myr, cooler (Ovi absorbing) and hotter (X-ray 
emitting) gas co- exist. This is a direct result of the turbu- 
lent and inhomogeneous density and temperature struc- 
ture in realistically evolved superbubbles. Turbulent mix- 
ing generates gas in the whole thermally unstable regime 
between 10^ - 10^ K. While the amount of Ovi in the 
LB has been measured fairly accurately (see Sect. I3.2|l . 
the source of the soft X-ray emission is more elusive. First 
of all, it is more difficult to localise, as there are most 
likely no dense clouds inside the LB, which can be used 
for shadowing experiments, thus allowing us to separate 
back- and foreground emission. Secondly, a yet uncertain 
amount of soft X-rays will be produced locally by solar 
wind charge exchange reactions with heliospheric plasma 
(Cravens 2000). And thirdly, the LB plasma need not be 
in CIE. It has been shown that if a plasma undergoes fast 
adiabatic expansion its kinetic temperature will drop, but 
recombination can not follow fast enough, thus mimick- 
ing an "overionised" plasma (similar to the solar wind), 
where recombination is delayed. In essence, the gas could 
be at a few times 10^ K, but still emit copiously in soft 
X-rays (Breitschwerdt & Schmutzlcr 1994). We are cur- 
rently working on these non-equilibrium ionisation (NEI) 
models, and will present results in a forthcoming paper. 

3.2. Ovi distribution in the solar neighbourhood 

We have determined the amount of Ovi along the lines of 
sight (LOS) for an ISM plasma in CIE with solar abun- 
dances (Anders & Grevesse, 1989; for details see AB05). 
To measure the Ovi column density distribution inside 
LB and LI, we took 91 LOS extending from the Sun and 
crossing LI (the hot pressured region 200 pc to the right 
of the LB) from an angle of —45° to +45° as marked in 
the top panel of Fig. ^ Fig. El (top) shows the history (for 
14.3 < t < 14.7 Myr) of the variation of the LOS averaged 
Ovi column densities, (N(Ovi)), with distance along these 
LOS. They decrease steeply with time after the last SN 
in the LB has exploded, from 4 x lO^'^ to 3 x 10^^ cm~^, 
because no further heating is taking place. For LOS sam- 
pling gas from outside the LB (i.e., Ilos > 100 pc) we 
have (N(Ovi)) > 2 x 10^^ cm^^. The histograms of col- 
umn densities obtained in the 91 LOS for t — 14.5 and 
14.6 Myr (Fig. El bottom) show that for t = 14.6 Myr all 
the LOS have column densities smaller than 10^^'^ cm~^, 
while for t = 14.5 Myr 67% of the lines have column den- 
sities smaller than lO^'^ cm~^ and about 50% of the lines 
have N(Ovi) < 7.9 x 10^^ cm~^. The maximum column 
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Fig. 2. Temperature distribution in the Galactic plane in the range 10'' < T < 10^ K of the LB and LI region, top: 
at t = 14.5 (left), t = 14.6 (middle) and t = 14.7 Myr (right), and bottom: at t = 14.8 (left), t = 14.9 (middle) and 
t — 15.0 Myr (right), after the first LB explosion. Note that in LI, SN explosions are still occurring (see bottom panel), 
while the LB is extinct. 



density is 1.78 x 10"'^'^ cm~^ in excellent agreement with the 
observations discussed in Savage & Lehner (2006), which 
is 1.1 times the average value of 1.6 x 10^^ cm~^ inferred by 
Shelton & Cox (1994) from a reanalysis of Copernicus 
absorption line data (cf. Jenkins 1978). 

In Fig. 0] we compare FUSE data by Oegerle et al. 
(2005; triangles) and Savage & Lehner (2006; circles) 
with simulated minimum (red squares), maximum (green 
squares) and averaged (blue squares) column density of 
Ovi measurements along the 91 LOS outlined in Fig. 1 
(top). It can be clearly seen in Fig. 0] that for t = 14.5 
Myr the calculated N(Ovi) distribution in the LB is sim- 
ilar to that observed with FUSE. These results allow us 
to constrain the age of the LB to be 14.5±q;4 Myr, and 
capture at the same time all the relevant LB properties, 
such as size, age, dynamical instabilities in the interaction 
shell, amongst others. 

Successive blast waves advance very fast through the 
low-density cavity, but slow down considerably later, as 
they run into the dense shell. As a consequence an asym- 
metric reverse shock leads to shear flow and turbulence 
inside the cavity, with a largest eddy size of a fraction 
of the bubble diameter, typically Z < 75 pc (Avillez & 
Breitschwerdt 2006). Thus the turnover time scale is about 
Tm ^ l/cs < 3.7 X 10^ yr, for an averaged sound speed 



of Cs « 200 km/s after a SN explosion. As the last SN 
occurred about 0.5 Myr ago, we expect that SN ejected 
oxygen inside the LB at present time should have a rather 
uniform O/H distribution, a result that is also confirmed 
by FUSE observations of the O/H ratio (Moos et al. 2002). 



4. Discussion and Conclusions 

According to Cox (2004) a crucial test for any viable LB 
model is the reproduction of the measured amount of Ovi. 
By embedding the LB into a typical inhomogeneous ISM, 
we obtain bubbles, which have more internal structure, 
exhibit on average lower pressure, consistent with that 
observed in the local clouds, show substantial mass load- 
ing (cf. Dyson & Hartquist 1987), and are more suscepti- 
ble to break-up of the surrounding shell (both in the disk 
and perpendicular to it), than was reported by previous 
authors. In fact there is some evidence from stellar ab- 
sorption line studies towards high latitude stars, that the 
LB might be a Local Chimney (Welsh et al. 1999). Break- 
out of gas flows from superbubbles might also be a clue 
to the solution of the "energy problem" observed in LMC 
bubbles (Oey & Garcia-Segura, 2004), which seem to be 
systematically too small compared to similarity solutions. 
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Fig. 3. Top: Averaged Ovi column densities (over 91 LOS 
between ±45° as marked in Fig. ^ top) as a function of 
LOS path length at 14.3 < t < 14.7 Myr of Local and LI 
bubbles evolution. Bottom: Histogram of the percentage 
of LOS as a function of N(Ovi) within the LB at f = 14.5 
and 14.6 Myr. 



Finally, we would like to stress that our LB evolution 
model has the highest resolution so far obtained, showing 
structures down to 1.25 pc. We have repeated our runs 
with 0.625 pc finest resolution to confirm that our results 
are not resolution dependent. This will be a prerequisite to 
model in a next step also the EUV and soft X-ray emis- 
sion, including the full non-equilibrium ionisation struc- 
ture. These calculations are underway. 
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